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We consider contact line deposition and pattern formation of a pinned evaporating thin drop. We 
identify and focus on the transport dynamics truncated by the maximal concentration, proposed by 
Dupont 1], as the single deposition mechanism. The truncated process, formalized as "pipe models", 
admits a characteristic moving shock front solution that has a robust functional form and depends 
only on local conditions. By applying the models, we solve the deposition process and describe 
the deposit density profile in different asymptotic regimes. In particular, near the contact line the 
density profile follows a scaling law that is proportional to the square root of the concentration ratio, 
and the maximal deposit density /thickness occurs at about 2/3 of the total drying time for uniform 
evaporation and 1/2 for diffusion-controlled evaporation. Away from the contact line, we for the 
first time identify the power-law decay of the deposit profile with respect to the radial distance. In 
comparison, our work is consistent with and extends previous results Q- We also predict features 
of the depinning process and multiple-ring patterns within Dupont model, and our predictions are 
consistent with empirical evidence. 



PACS numbers: 68.15.+e, 47.20.-k, 68.60.D\ 



I. INTRODUCTION 



Evaporative contact line deposition, also known as the 
"coffee-drop effect" has been the subject of several recent 
papers [2], |3|, |4|, la la [3] ■ The physical problem originates 
from a simple phenomenon of everyday life: when a drop 
containing a solute such as coffee dries on a surface, the 
solute is driven to the contact line and forms a charac- 
teristic ring pattern. The basic mechanism behind the 
phenomenon is now well understood : the drying drop 
is pinned at the contact line so that there must be an 
outward capillary flow to replace the liquid evaporating 
at the edge, and this flow carries solute particles to the 
contact line where the deposit is formed. 

This simple phenomenon is important in many scien- 
tific and industrial applications [1, H, Ha EH, E2, ED, EH . 
The evaporation mechanism can create very fine lines of 
deposit in a robust way that requires no explicit form- 
ing. It can be used to concentrate material strongly in a 
controllable way. It also creates capillary flow patterns 
which are useful for processing polyatomic solutes like 

dna [H,E1,E3- 

The deposition process and the resulting patterns are 
affected by various physical and chemical conditions. The 
wetting property of the substrate and surface roughness 
affect the geometry of an evaporating drop dUEll]- Also, 
within the drop there can be a convective flow induced 
by gradient in temperature, surfactant concentration, or 
surface tension [2(J, HH [H, H3 1 and it can lead to very 
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complicated, even fractal-like, deposit patterns. Further- 
more, some mechanisms, such as contact line depinning 
[E 0], are s tiH n °t clear, and there are no satisfactory 
explanations for multiple-ring deposit patterns [25l |26| . 

Several models have been established under the com- 
mon assumptions of slow evaporation and creeping flow 
so that quasi-equilibrium dynamics is warranted: the dy- 
namics is determined by the time evolution of the equilib- 
rium properties of the evaporating system. These models 
use different conditions and assumptions on such proper- 
ties as geometry, evaporation profile, and deposition cri- 
teria. In addition to circular thin drop, there have been 
studies on contact line deposition in an angular region 
[27l [H, , where the nontrivial geometry at the tip of 
the angle induces a crossover in deposit properties. Both 
uniform evaporation flux and diffusion-controlled flux, 
which diverges at the contact line, are considered, and 
they yield different deposit density profiles and growth 
properties along the contact line in this asymmetric ge- 
ometry. 

In a circular-symmetric geometry, such as a circular 
evaporating drop, the deposit along the contact line is 
always uniform. It is thus more important to understand 
in this case the spatial deposit profile and growth dy- 
namics in a dimension perpendicular to the contact line, 
that is, how the solute particles form the deposit and ac- 
cumulate toward the center of the drop from the contact 
line. Dupont [l[ suggested there be a maximal volume 
concentration C max for solute particles such that once 
the concentration increases to C max , the horizontal trans- 
port of solute particles stops, and those particles form the 
deposit. Following Dupont's deposition criterion, Popov 
01 studied a model with a pinned round evaporating thin 
drop and diffusion-controlled evaporation flux. There are 
two regions in his model: a deposition region, where the 
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FIG. 2: Direct numerical simulation of the deposit pattern. 
The areal density profile is represented by variation in dark- 
ness. (Courtesy Todd Dupont) 



FIG. 1: Sketches of models of a pinned evaporating drop. Fig- 
ure A shows Popov's model, where the deposit (shaded area) 
maintains its shape and modifies the geometry of the remain- 
ing evaporating drop. Figure B shows the model considered 
in this work. The deposit profile does not affect but passively 
changes with the geometry of the evaporating drop. Figure 
C shows the projected picture from above. The deposit (dark 
area) grows from the contact line toward the center of the 
drop. 



solute concentration is C max ; a transport region, where 
the solute concentration is smaller than C max and hori- 
zontal solute transport is allowed. Popov assumed that 
the deposition near the contact line affects the geome- 
try of the drop and the evaporation-induced flow field 
(Fig 1.2A). As a result, Popov's model is characterized 
by several variables described by a system of coupled dif- 
ferential equations, of which a complete solution is not 
available. 

In Popov's picture, the deposit thickness increases 
monotonically toward the center of the drop until an 
abrupt vertical wall is formed at the end of the total 
drying time. The scale of the thickness is proportional 
to -\/^0i where So is the concentration ratio defined as 
the initial uniform volume concentration Co divided by 
C max . This \/So-dependence near the contact line had 
been known before [f| and is understandable heuristi- 
cally. The volume of the deposit, which is proportional 
to the square of the thickness, multiplied by the con- 
centration C max , must give the total solute mass in the 
deposit that is proportional to Co- However, as this de- 
pendence is largely asymptotic toward the contact line, 
where the time dependence is suppressed and the hy- 
drodynamics and the deposition dynamics are largely 
separated, Popov's coupled system may not be neces- 
sary to derive this result. Moreover, the complexity in 
Popov's mathematical formulations may inadvertently 
obscure the simple and fundamental underlying physi- 
cal process. Moreover, a realistic deposit clearly follows 
a rich profile that shows different scaling properties in 
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FIG. 3: Areal density profile E with Co/C ma x = 0.002, 
where Co = 0.001 is the initial uniform volume concentra- 
tion throughout the drop. The dotted line is the initial areal 
density identified as the initial drop height profile multiplied 
by Co- The solid line is the final areal desntity (Courtesy 
Todd Dupont) 



different regions. Thus to interpret the deposit proper- 
ties within a single asymptotic regime and with a unique 
yJo is n °t satisfactory. 

Following the truncated transport dynamics and the 
simpler assumptions that the hydrodynamics and the de- 
position process are fully separable, Dupont lj showed 
by direct numerical simulation a ring-like deposit pattern 
with variance in areal density along the radial direction 
(Fig. [2] and Fig. [3]). The simple mechanism thus pro- 
duces a realistic deposit pattern with a pronounced maxi- 
mal density near the contact line. In particular, the areal 
density shows a potentially analyzable characteristic pro- 
file that manifests multiple asymptotic regimes, and that 
could reveal robust properties of contact line deposition 
under general conditions. 

In this work we re-visit Dupont 's criterion with mini- 
mal restrictions and simpler assumptions. We consider a 
pinned circular thin drop. In particular, we assume that 
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the deposition process does not interfere with the hydro- 
dynamics such as the drop geometry and flow velocity 
field. We consider both the uniform evaporation and the 
diffusion-controlled evaporation conditions, to be defined 
below. To interpret the deposit density, we assume that 
the solute particles in the deposition region, where hor- 
izontal transport is stopped, can still move vertically so 
that the thickness of the deposit is determined by the 
geometry of the drop and changes with time. The fi- 
nal deposit profile is thus described by its areal density 
when the drop dries up and the thickness becomes zero 
throughout the drop (Fig. 1.2B). To interpret the final 
deposit in terms of the areal density profile is not essen- 
tial to the model, however, as we shall discuss later that 
the model is also valid for a deposit with finite thick- 
ness as long as the deposition does not interfere with the 
hydrodynamics of the evaporating drop. 

We thus want to follow a reductive approach as we be- 
lieve that the fundamental and robust physical properties 
underlying complex phenomena should be captured by a 
simple mathematical structure with minimal restrictions. 
By this approach, we want to achieve two goals. First, 
instead of being exhaustive in considering all the pos- 
sible conditions, we want to understand the robustness 
of the model: among the various factors that character- 
ize the problem what are most important in determining 
the key deposition properties and in yielding physically 
feasible and realistic deposit patterns. Second, what is 
the appropriate mathematical structure to describe those 
dominant factors, and can this mathematical structure be 
generalized to other physical systems. 

In the rest of the work we first establish in Section 2 a 
one-dimensional simplification of the deposition process 
that we denote as "pipe models" . These models capture 
the important characteristics of the truncated transport 
dynamics. We then review the major results of the hy- 
drodynamic properties of an evaporating drop and apply 
the pipe models to the evaporative deposition process 
with uniform evaporation in Section 3. In Section 4, we 
solve for the areal density profile and identify the scal- 
ing laws in different asymptotic regimes. We consider 
in Section 5 the deposition properties with the diffusion- 
controlled evaporation. In Section 6, we compare our 
findings with Popov's, extending our results to interpret- 
ing the deposit in terms of the thickness profile. We 
then make a few comments on depinning and formation 
of multiple-ring patterns within the truncated transport 
dynamics scheme in Section 7. We discuss experimental 
issues and possible generalizations of the model in Section 
8, and present the conclusion in Section 9. In addition, 
we provide an alternative derivation of the equation of 
motion for the shock front in Appendix. 



II. PIPE MODELS 

The transport process, truncated by the maximal con- 
centration C max , is naturally divided into two regions: a 




FIG. 4: A uniform pipe. Solute particles with uniform volume 
concentration Co is carried at a constant velocity v toward the 
left end of the pipe, where deposit is formed with concentra- 
tion C m ax > Co- l(t) is the position of the moving boundary 
that separates the two regions with different volume concen- 
tration. 



transport region where the concentration is smaller than 
Cmax: an d a deposition region where there is no horizon- 
tal transport of solute particles. Solute mass is carried 
from the transport region to the deposition region, and 
the growth of the deposit is described by a characteristic 
moving front that separates the two regions. 

We establish in this Section the so-called "pipe mod- 
els" to formalize the truncated dynamics under different 
conditions. Although some predictions made by these 
models depend on pipe specifications, we find that certain 
important properties, such as the shock front velocity, are 
actually model-independent and thus generalizable. We 
shall show later how these universal properties may help 
explain the observed robustness in the general deposition 
phenomena. 

To simplify mathematical formulations, in the follow- 
ing models we only treat one-dimensional case, or simi- 
larly radial symmetric case. However, we believe the key 
findings are readily applicable to higher dimensions. 



A. A Uniform Pipe 

We consider a uniform pipe with simple properties 
(Fig. 2]): incoming solute particles with constant vol- 
ume concentration Co are transported at a constant ve- 
locity v to the end of the pipe where the carrier fluid 
exits the pipe but the solute particles remain to form a 
deposit. The solute particles stop moving when their con- 
centration exceeds C max , which therefore is the volume 
concentration of the deposit. To focus on the horizontal 
transport properties, we ignore the variance in velocity 
along the cross section of the pipe. This simplification 
is compatible with the basic assumptions made for the 
evaporative deposition problem as we shall show in the 
next Section. 

The configuration is thus divided into two regions: a 
transport region with concentration Co and a deposition 
region with concentration C max . The boundary that sep- 
arates the two regions is characterized by its position 
l{t). Local conservation of solute mass at the boundary 
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demands: 



dl 
dt 



Co 



C n 



Co 



1 



(1) 



where S = C /C max . 

Several observations immediately follow from Eq. (fTj). 
First, the boundary moves at a velocity proportional to 
the transport velocity v, and the proportional coefficient 
is determined by the ratio of the two volumes concen- 
trations. Second, there must be a finite gap in volume 
concentration across the boundary, i.e., S < 1, other- 
wise the velocity dl/dt would diverge, which is not physi- 
cal. The moving boundary is therefore a shock front that 
represents a discontinuity in local volume concentration 
profile. 

It is worth noting that although the two-component 
simplified geometry with constant Co and C max in the 
first place warrants the discontinuity and leads to the 
shock front, the implications of Eq. Jl} are more general. 
If the transport process is characterized by the trunca- 
tion criterion and the two-component configuration and if 
the boundary that separates the two regions moves with 
a finite velocity, there gap in concentrations across the 
boundary must exist as we now show. 
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FIG. 5: A realistic pipe, r is the radial distance. The height 
profile of the pipe is h(r,t). 2nrh(r,t) characterizes the cross 
section area, which is not shown, to account for the circularly 
symmetric geometry. The concentration profile C(r, i) in the 
transport region, the flow velocity field v(r, t), and the critical 
concentration C max (r, t) are allowed to change with both the 
radial distance and time. l(t) is the boundary that separates 
the deposition region and the transport region. In a small 
time interval At, l(t) moves to l(t + At) — 1(f) + Al. Am is 
the change of solute mass in the transitional region. 



B. A More Realistic Pipe: Non-Singular Case 

A typical evaporating drop that we consider has a cir- 
cularly symmetric geometry and an outward flow field 
along the radial direction. To adapt to this specific ge- 
ometry, in a more realistic pipe model (Fig. [5J we in- 
troduce h(r, t) as the pipe's height profile, where r is the 
radial distance, and allow other characteristic quantities 
such as the maximal volume concentration C max and the 
transport velocity v to vary both in r and time t. In ad- 
dtion, to account for the circularly symmetric geometry, 
our pipe must have a cross section area that varies along 
the radial direction as 2irrh(r,t). 

This manipulation with respect to the circularly sym- 
metric geometry may seem greatly constraint the pipe 
models. However, we shall show below that the major 
predictions made from the pipe models are actually in- 
dependent of the specific form of the cross section area 
as long as it does not introduce significant singularities. 
The term 2irrh(r, t) has been introduced only to conform 
with the drop geometry which is not really a "pipe" . 

We treat in the real pipe h, C max , and v as independent 
specified physical quantities that have independent dy- 
namics. This is not necessarily true in real settings such 
as the evaporative deposition we shall consider later, but 
it will allow us to derive general results and will pose not 
difficulties to implementing these results to specific real 
situations where the dynamics may couple. There are 
perhaps some conceptual difficulties with the functional 
form C max (r, t). For example, because it evolves in time, 
what if for a position in the deposition region the C max 
becomes smaller than the actual volume concentration 



in a later time? We do not consider this special case 
as it does not occur in the evaporative deposition prob- 
lem. Nonetheless, it can be relevant in other systems 
with more complicated transport properties and is cer- 
tainly an interesting problem for its own sake. 

We again impose the truncation condition that the hor- 
izontal transport of solute mass is truncated once the lo- 
cal volume concentration C(r,t) reaches C max (r,t). As 
in the uniform pipe model, we consider the position l(t) 
of the moving front that separates the deposition region 
and the transport region. To simplify mathematical anal- 
ysis, we assume that all the introduced functions are at 
lease twice diffcrcntiable and that there are no singular- 
ities. We shall however discuss effects of discontinuities 
and other singularities later. 

In the deposition region there is no horizontal trans- 
port of solute mass, and local conservation of mass re- 
quires C(r, t)h(r, t) to be time-independent. In the trans- 
port region, solute mass satisfies the equation of continu- 
ity 

- (Chr) + ^ (vChr) = 0. (2) 

As shown in Fig. during the time interval At the 
shock front moves from l(t) to l(t + At) = l(t) + Al. 
To be consistent with the later-introduced context of an 
evaporating drop, we choose Al < 0. l(t) thus decreases 
in time and moves toward the right end of the pipe. 

We consider the change of solute mass Am in the tran- 
sitional region [l(t) + Al, l(t)]: 
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Am = Ami 

rt+At 



2itI{t) 



x [C max (l (r) , r)h(l (r) , r) - C{1 (r) , t)/i(Z (r) , t)] 

x (-§)*, (3) 

i.e., Am is equal to the amount of solute mass in the 
region when it just becomes a part of the deposition re- 
gion less the amount of mass in it when it is still in the 
transport region. 

We expand Ami up to the second order of At: 



Ami = -2nl(t) 



dl 



x [C max (l(t),t)h(l(t),t) - C{l{t),t){h{l{t),t)] -At 

dt 

1 d ,„ . dl 
x [C max (Z (r) , r)h(l (r) , r) - C{1 (r) , t)h{l (r) , *)] }| T=t 



where 5 = 5{r, t) = C{r, i)/C max (r, t). 

Eq. has the same form as Eq. |T]) in the case 
of the uniform pipe though the constant quantities have 
been replaced here by local values of general functions 
dependent on position and time. The moving front is a 
shock front with 5 < 1. The front velocity is proportional 
to the local transport velocity with opposite direction, 
and the proportional coefficient is determined by the local 
critical concentration C max and the local concentration C 
in the transport region. Eq. is thus the general non- 
singular form of the equation of motion for the shock 
front in the truncated dynamics. 

It also shows that the discontinuity condition 5 < 1 is 
not due to the instantiations of constant Co and C max but 
the requirement of finite velocity and is inherent in the 
truncated transport dynamics. As both concentrations 
are given functions, Eq. (jSJ) imposes a specific moving 
boundary l(t) such that a time-dependent jump in con- 
centration across the boundary is properly maintained 
and the local solute mass is thus conserved. 



x (A<) 2 + 0{{Atf) 
= -2vr/^(C ma 



C^At 
' dt 



dt 

0((Atf) 



C)hl 



{Aty 



C. A More Realistic Pipe: Singular Case 

If the local height profile h(l(t), t) is zero, the first or- 
der term in At vanishes for both Ami and Am-2, and the 
second order terms in At thus must match. The condi- 
tion Am = Ami = Am,2 requires: 



(4) 



Alternatively, Am is also equal to the amount of solute 
mass transported through l(t + At) — l(t) + Al from time 
t to t + At: 



Am = Am.2 

r t+At 



-i; 



2irl(t + At)C(l(t + At),r)h (l(t + At), r) 



and 

dl 
dt 



d ( ' „ N , ,dl\ d , .,dl 



2 (d/dl) (IChv) (dl/dt) + (d/dt) (IChv) 



(d/dt) {{C n 



C) hi) + I {d/dt) {Ch) 



(9) 



(10) 



x v(l{t + At),r)dr. (5) 
We also expand Am-2 up to the second order of At: 



Am 2 = 2nlChvAt + 2ir 
x {At) 2 + 0{{Atf). 



(6) 



The general form Eq. (|10p can be simplified under 
certain physical conditions. In particular, we consider 
the situation where h{l{t),t) = dh{l{t),t)/dt = and 
dh/dl 7^ 0, which corresponds to the case when the shock 
front reaches a bottleneck of the pipe. In addition we 
assume all other quantities are nonvanishing at the bot- 
tleneck. Thus at l{t) 



We can now solve for dl / dt given the condition Am = 
Ami = Am 2 . If h{l{t),t) is nonvanishing, the terms in 
the first order of At must match, and thus: 



2%lh{C max ~C)j t = ^nlChv, (7) 



m _ c{i{t),t) s_ 

dt~~ v{ C max {l{t),t)-C{l{t),t) - ~ V T~5' 



-{IChv) = Cv-{hl), 



at {Ch) = °' 



(ii) 



(12) 



and 



^ ((C max - C) hi) = 1 ((C max - C) hi) 
dt ot 



(8) 



+ ^ C max - C) hi) — = C n 

ol dt 
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Also, since C(l(t),t) satisfies the equation of continuity 
(Eq. ©), 



(d/dt) (IChv) = -v 2 C^ (hi) 



(14) 



Substituting Eqs. (fTTj) to (fl"4")) , the shock front velocity 
Eq. (fTO)) is reduced to 



dl 
dt 



C 



C n 



c 



dl/dt 



It can be solved for dl/dt: 

dl _ 
di ~ 



V6 



i-Vs' 



(15) 



(16) 



where S = 5(l(t),t) = C(l(t),t)/C max (l(t),t) is the local 
concentration ratio at l(t) such that h(l(t),t) vanishes. 

The functional form of Eq. (fl6|) is similar to Eq. ([_]) 
for the nonsingular case with 5 replaced by \fS. How- 
ever, it is worth noting that Eq. (TT6")) holds only at those 
discrete points along the pipe where the height is zero. 
In this sense, Eq. (|16p is not a differential equation but 
a boundary condition or a series of boundary conditions 
determined by the properties of the pipe in question. We 
shall come back to this point in Discussion. 



D. Truncated Transport Dynamics and the Shock 
Front 

We have shown that the truncated dynamics in various 
pipe models is described by a moving shock front l(t), of 
which the equation of motion is determined by the local 
values of the characteristic quantities at the front under 
the condition of mass conservation. We now establish 
the formal mathematical framework for one-dimensional 
truncated transport dynamics and derive the equation 
of motion for the characteristic shock front directly. We 
only treat non-singular case and omit rigorous discussion 
of boundary and initial conditions. In particular, we do 
not consider the problem how the dynamics is initiated 
but once it has been established what the equation of mo- 
tion should be. The solution to the former depends on the 
boundary and initial conditions but is less relevant to ap- 
plying the models to the evaporative deposition problem. 

We consider a one-dimensional transport equation 
which imposes both the continuity condition with respect 
to the generalized mass ip(r,t)C(r,t) and the truncation 
condition imposed on C (r, t) : 



0_ 

dt 



(<p(r,t)C(r,t)) 

d 



+ — (v(r,t)<p(r,t)C(r,t)6(C max (r,t) - C{r,t)) = 0, 

(17) 

where v(r, t) is the velocity field, C max (r, t) is the trun- 
cation threshold, and 9 is the step function defined by: 



c 






QJr,t) 








| ________ 



lit) 



r = 



FIG. 6: Sketch of the truncated dynamics solution at time 
t. The solution consists two parts: C(r, t) in the transport 
region and C(r,i) in the deposition region. The two regions 
are separated by the shock front l(t). There is a nonzero jump 
in C(r,t) across the shock front. 



6(x) = i 1 ' X> °> 
y ' \0, x < 0. 



(18) 



For general purposes, tp(r,t), C (r, t), and C max (r, t) 
in Eq. ([17]) have no specific physical meaning attached 
though (p(r,t) can be regarded as the ratio between the 
generalized mass and the quantity on which the trunca- 
tion condition is imposed. In the case of the pipe models 
we have considered, C(r, t) is the solute volume concen- 
tration and ip(r,t) = rh(r,t) where h(r,t) is the height 
profile of the pipe. For non-singular solutions, we as- 
sume (p(r,t), v(r,t), and C max (r, t) are positive and at 
least twice differentiable. 

We seek solution that admits a moving shock front l[t): 

C(r, t) = C(r, t) + 0(l(t) - r)(C(r, t) - C(r, t)), (19) 

where C (r, t) and C (r, t) are again at least twice differen- 
tiable, C(r,t) < C max (r,i) for r > l(t), and C(l(t),t) = 
C m ax- We show the sketch of C(r, t) in Fig. _] 

Substituting Eq. (TT9")) into Eq. (fT?)) . all singular terms 
must cancel each other. Thus, 



(y{C-C)j t +wC)5{l-r) 

and, given (p(l(t),t) is not zero, 

C(l(t),t) 



0. 



dl 
di 



-w(r, t) 



C max (l(t),t) - C(l(t),t) 



(20) 



(21) 



Eq. (f2"Tj) is the general result of pipe models, and in the 
non-singular case it does not depend on tp(r,t). 

The coffee stain problem we shall discuss in the next 
Section corresponds to the specific case where ip(r, t) = 
rh(r,t), C max (r,t) = constant, and C(r, 0) = Cq. 
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III. EVAPORATIVE DEPOSITION WITH 
UNIFORM EVAPORATION 



The evaporative deposition process is a special case of 
the pipe models above: hydrodynamics of the thin evap- 
orating drop determines a height profile h(r,t) and gen- 
erates a flow velocity field v(r,t); solute particles with 
initial uniform volume concentration Co are carried by 
the evaporation-generated flow toward the contact line, 
and the transport process is truncated by a maximal al- 
lowed concentration C max ; once the volume concentra- 
tion reaches C max , the horizontal transport is stopped 
and deposit forms. 

The volume concentration profile of solute particles is 
thus divided into two regions: 

Transport region. Solute particles move with the 
evaporation-induced flow at the same velocity and are 
transported toward the contact line along the radial di- 
rection. The volume concentration in this region solves 
the continuity equation with a finite velocity field. 

Deposition region. The horizontal transport of solute 
particles is truncated, and only vertical movement is al- 
lowed afterward. The solute mass at each position, i.e., 
C(r,t)h(r,t) is thus conserved in time. All the solute 
particles eventually deposit at the position and form the 
areal density profile at the total drying time. We also as- 
sume that the deposition process is decoupled from the 
hydrodynamics of the evaporating drop, and both the 
height profile h(r,t) and the velocity field v(r,t) are not 
affected by the deposition process. 

There must be a moving shock front l(t) that separates 
these two regions. The velocity of the front dl/dt must 
obey the predictions made by the pipe models. 

To apply the pipe models, we need to specify the height 
profile h(r,t), the velocity field v(r,t), and the solute 
volume concentration in the transport region. Most of 
these have already been explored in the existent litera- 
ture |1, 0] though some exact results are not immediately 
available. For the sake of completeness, we shall nonethe- 
less go through the whole derivation process from our own 
perspective. 




FIG. 7: Sketch of an evaporating thin drop. We assume the 
drop is pinned during the time under consideration. The ra- 
dius of the drop R is thus a constant. The geometry is de- 
scribed by a height profile h(r, t) that depends on the radial 
distance and time. 8(t) is the contact angle that decreases lin- 
early in time in our approximation. J(r) is the evaporation 
flux profile. v(r,t) is the flow velocity field that also depends 
on the radial distance and time. 



pressure difference across the liquid-air interface. 1 The 
shape of the drop is thus determined by 



Ap = -2Ha, 



(22) 



where Ap — Ap{t) is the pressure difference, H = H(t) 
is the mean curvature of the drop surface, and cr is the 
surface tension, ft is well known that for a circular drop 
the geometry that solves Eq. (|2"2"|) is a spherical cap with 
the contact angle 6 determined by: 



sin0(t) = 



R Ap{t) 



2 a 

The height profile h(r,t) in this case is 



(23) 



h(r,t) 



y sin 2 6{t) 
R 2 -r 2 
2R 



r 2 - R cot 9{t) 
6(t), (9<1. (24) 



If the evaporation flux does not change with time, it is 
confirmed experimentally [f| that when is small, 2 the 
contact angle 9(t) decreases linearly with time for the 
most part of the drying process: 



Review of Evaporating Thin Drop 
Hydrodynamics 



We consider a circular evaporating thin liquid drop 
(Fig. [7]). The drop is pinned in the drying process, and 
the drop radius R is thus a constant. The geometry of 
the drop is determined by the contact angle 0(t) and the 
height profile h(r, t). 

Under the assumption of slow evaporation, the quasi- 
equilibrium condition can be applied so that at any time 
the drop maintains mechanical equilibrium against the 



9(t) = 9i(l--), 



(25) 



where 8i is the initial contact angle and tf is the total 
drying time that depends on specific forms of evaporation 



1 This condition hoids as long as the flow velocity within the drop 
is much smaller than the characteristic velocity u/3r] where a is 
the surface tension and v is the dynamic viscosity 7]. For water 
the characteristic velocity is about 24 m/s. 

2 The initial contact angle 0; in Deegan's experiments was about 
0.1 to 0.3 radians 
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profile. In what follows we shall restrict our discussion 
to thin drops so that the linear relation Eq. (f^5|) is war- 
ranted. 

Loss of solvent due to evaporation generates inside 
the drop a flow field u(r, z, t) with horizontal component 
U a [r, z, t), where z represents the coordinate along the 
vertical direction. We introduce a vertically averaged ve- 
locity field 



v(r,t) 



1 



h(r,t) 



u s (r, z, t)dz. 



Hr,t) jo 

Local conservation of solvent mass then demands 
J(r) 



(26) 



V • (hv) 



1 + (Vhy + d t h = 0, (27) 




FIG. 8: Sketch of the definition of m(r,t). A solute particle 
at ro at time zero moves to r at time t. ro as a function of r 
and t is given by Eq. (|33[1 . m(r, t) is equal to the total solute 
mass between ro and r (the shaded area) at time zero. 



where p is the density of the solvent, and J(r) is the 
evaporation profile defined as mass loss per unit projected 
area per unit time. For a thin circular drop v(r, t) is along 
the radial direction and \Vh\ <C 1, and Eq. ([27]) can be 
integrated to solve for v: 



v(r,t) 



1 

rh ./,, 



J 



+ d t h ) rdi 



(28) 



For uniform evaporation in which J(r) = const., it can 
be shown that: 



/ % lr / t 



(29) 



where the total drying time t / , when the whole drop dries 
up via evaporation, is 



ttR 2 J 



2irh(r, 0)rdr 



e l R P 

4J ■ 



(30) 



To simplify mathematical formulations, we choose R as 
the unit of length, tf as the unit of time, and R/tf as the 
unit of velocity. Thus in this model the dynamics of an 
evaporating drop is characterized by two dimensionless 
functions: 



h(r,t) = j(l-r 2 )(l-t), 



(31) 



B. Solute Volume Concentration Profile 

1. Concentration Profile in the Transport Region 

We now derive the volume concentration C (r, t) in the 
transport region, assuming the defined height profile (Eq. 
([31])) and velocity field (Eq. ([32]) ). 

For radial distance r and time t we define ro (r, t) to 
be the initial position such that solute particles start at 
time zero from ro reach r at time t. tq can be solved via 
integrating the velocity field (Eq. (|32l ): 



r (r,t) =r(l-t) 



(33) 



We next define m(r, t) to be the amount of solute pass- 
ing through r from time zero to time t (Fig. [5]). Then 
m(r, t) is equal to the total amount of solute between r 
and ro(r, t) at time zero: 



m{r,t) = 2ir C h(r' ,0)r'dr 

Jr (r,t) 
2 

tX\a L.(i-(i-t) 



^ - U 



(34) 



We consider a ring-like region between r and r+ Ar. At 
time zero the total amount of solute mass in this region 
is 



v (r,t) = -r(l - t)' 



(32) 



where we have used the same notations r and t to repre- 
sent the dimensionless distance and time respectively for 
simplicity. 

In deriving v(r, t) (Eq. (|3"2"|) i, we have used the condi- 
tion that the mass of the solvent is conserved locally (Eq. 
([27]) ). and this condition involves h(r, t). Thus v(r, t) and 
h(r, t) have coupled dynamics as opposed to the general 
situation we considered in the real pipe models. This 
coupling will not pose any difficulty, however. 



m = 27rrAr/i(r,0)C . 
At time t, it becomes 

rat = 2irrArh(r,t)C(r,t). 



(35) 



(36) 



The change of solute mass up to time t in this region 
should be equal to the total amount of mass moving into 
r less the total amount of mass moving out of r + Ar, 



Ar« = m(r, t) — m(r + Ar, t). 



(37) 
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i i , 1 , 1 , 1 , 1 , 1 

0.0 0.2 0.4 0.6 0.8 1.0 



FIG. 9: The solute volume concentration profile C(r,t)/Co 
(Eq. (|38|ll in the transport region at t = 0.2, 0.4, 0.6, and 
0.8. 



that is, the change of the total mass is equal to the 
amount of mass injected into the drop (at r = 0) less 
the amount of mass flowing out of the drop boundary 
(at r = 1) that immediately forms deposit at the contact 
line. The transport process thus must be truncated with 
some finite threshold C max so that a moving boundary 
is formed and the deposit grows from the contact line 
toward the center of the drop. 

2. Concentration Profile in the Deposition Region 

We are now ready to apply the pipe models. There 
is a shock front l(t) that separates the transport region 
and the deposition region. The equation of motion for 
the shock front is given by Eq. for the nonsingular 
case, where 6 — C(r,t)/C max is the concentration ratio 
at the front. Substituting the height profile h(r,i) (Eq. 
(|3ip ). the flow velocity field (Eq. (|32|) ). and the volume 
concentration profile in the transport region (Eq. (|38|) ). 
the velocity dl/dt is given by: 



Conservation of solute mass demands mo 
thus 



Am 



mt, 



C(r,t) = C(r,t) 



lim 

Ar— >0 



Co 



h(r, 0) m(r,t) — m{r + Ar, i) 



= C 



(1-t) 



h(r,t) 
-V 2 _ r 2 



2nrArh(r, t) 



1 



(38) 



where, in accordance with the notations in pipe models, 
we use C (r, t) to denote the volume concentration in the 
transport region. As shown in Fig. [5J C(r,t) increases 
monotonically with both r and t. 

The concentration profile (Eq. {38])) is a "leaking so- 
lution" in the sense that the total amount of solute mass 
in the transport region is not conserved, and rightly so. 
This can be demonstrated in an ideal evaporative depo- 
sition process such that all the solute particles are even- 
tually carried to the contact line without truncation in 
volume concentration i.e., in the limit C max — > oo. It 
is verifiable that C(r, t) satisfies the continuity equation 
(Eq. ©): 

The total amount of solute mass at time t in the body of 
the drop is 

M(t) = [ 2-KChrdr. (40) 



rZr 



Thus the rate of change of the total solute mass is, 

dM(t) f 1 d , ~ s f 1 d 

1 2w—(Chr)dr = -2ir 



dt 



(i 



Or 



(yChr) 



= 2-KvChi 



r=0 



(41) 



dl 
dt 



I 



4(1 -t) gm~(l-P) . _l' 



(42) 



To determine l(t) completely, in addition to Eq. 
we also need to specify initial conditions. At time zero, 
the shock front l(t) starts at the contact line: 1(0) = 1. 
Also, physically dl/dt must have a well-define value at 
time zero, but naively for t = and 1(0) = 1 the right 
side of Eq. (j42"]) is not well defined. The pipe models 
again can shed some light on this. At contact line the 
height is zero, and the pipe models in the singular case 
(Eq. (Uni)) require 



dl 
~dt 



-u(0,0)- 



1 



1(0) 



1. 



t=o 



1-VSo 4l-V3o 

(43) 

where 6q — Co/C max . Eqs. (|42)) and (|43f thus completely 
define the time evolution of the shock front l(t). 

The seeming discontinuous "jump" in dl/dt from the 
V^-dcpendence at the contact line to (^-dependence once 
I moves away is artificial. The concentration profile in the 
transport region C(r,t) (Eq. I[38p) has no definite limit 
at the contact line for r — 1 and t — 0. Specific limit thus 
must be calculated along a definite path r — l(t). The 
transition from \/~5$ to S is actually continuous along this 
path. This can be demonstrated directly (see Appendix 
for details), and will be further clarified later when we 
discuss the deposit profile properties near the contact 
line. 

Once l(t) is known, the concentration profile in the 
deposition region can be determined under the condition 
that C(r,t)h(r,t) is conserved in this region. Thus for 
r > l(t) 



C(r,t) = C(r,t) = C n 



h(r,l-\r)) 
' h(r,t) 



(44) 
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where Z _1 is the inverse function of I: l _1 (r) gives the 
time when the shock front arrives at the position r. 

Combining the concentration profile C (Eq. (|38|) ) in 
the transport region and the concentration profile C (Eq. 
(|44p ) in the deposition region via the moving the shock 
front l(t), the total concentration profile is 

C(r, t) = C(r, t) + 9(r - l(t))(C(r, t) ~ C(r, t)). (45) 

Alternatively the shock front velocity (Eq. (|4"2|) ) and 



the boundary condition (Eq. ([43]) ) can be derived di- 
rectly via solute mass conservation. We show the deriva- 
tion in details in Appendix. 



c 



\ 

\ 

\ 

\ 

\ 
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\ 

\ 






\ 

\ 





r = l l(t) l(t) 



r = 



C. The Shock Front 

The region boundary l(t) is a moving shock front, and 
there is a finite gap in volume concentration across the 
boundary. 

Naively we can introduce a virtual moving front l(t) by 
extrapolating the concentration profile (Eq. ([55]) ) in the 
transport region with the condition C(l(t),t) = C max . 
Thus 

1/2 



l(t) 



1 



*0 



1 - S 1-S 



(l-t) 



-1/2 



(46) 



where <5 = Co/C max . l(t) represents the location of the 
moving front if there were no discontinuity in volume 
concentration. Evidently t < 1 — Sq for Eq. (PlU)) to be 
meaningful, that is, the virtual front l(t) arrives at the 
center of the drop before the whole drop dries up at t = 1 . 
We can also derive the equation of motion for l(t): 



dl 



= -v(l(t),t) 



lit 



(l-t) 1/2 -S 



1 So 



4 1 -So 



(47) 



(48) 



Thus the virtual front moves slower than the shock front 
starting from the contact line. 

Since for < t < 1 the volume concentration C(r, t) 
(Eq. (|38p ) in the transport region increases with r and 
at the shock front C(l{t),t) < C max = C(l(t),t), thus 

l(t) < l(t), (49) 

for < t < 1 - Si. 

We know that 1(0) = 1(0) = 
later that 1(1 - S%) = 1(1 - S%) 
and the virtual front thus start out from the contact line 
at time t — with the shock front moving faster. The 
virtual front then catches up, and both fronts arrive at 
the center of the drop at time t = 1 — 5q. The sketches 
of these two fronts at time t are shown in Fig. [101 an d 
their time evolutions are shown in Fig. [TTJ We note that 
both fronts spend most of the total drying time near the 
contact line (I — > 1). 



1, and we shall prove 
= 0. The shock front 



FIG. 10: Sketches of the shock front and the virtual front at 
time t. The shock front is ahead of the virtual front: l(t) < 
l(t). There is a nonzero jump in volume concentration across 
the shock front: C max — C(l(t),i), with C max = C(l(i),t). 




FIG. 11: Time evolution of the two fronts with <5o = 0.002. 
l(t) and l(t) start from the contact line (r = 1) at time zero 
and arrive at the center of the drop (r = 0) at time t = 1 — 6q . 



IV. DEPOSIT DENSITY PROFILE WITH 
UNIFORM EVAPORATION 

A. Areal Density 

When a drop dries up, a ring-like stain is left with 
nonzero width and variation in density. The spatial pro- 
file of the deposit pattern is characterized by the areal 
density S, defined as the amount of solute mass per unit 
area. For a pinned circular drop without irregularities, 
£ is a function of the radial distance r. 

Since in the deposition region there is no horizontal so- 
lute transport, at any position in this region all the solute 
particles distributing vertically above the position along 
the height of the drop will eventually deposit and form 
the areal density at the position. We consider a position 
with radial distance r from the center of the drop. At 
time l _1 (r) when the shock front arrives at r, the drop 
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2Z(r) 




0.2 0.4 0.6 0.8 1 



FIG. 12: The areal density profile 2£(r)/(0 i C Ina3t ) from the 
center of the drop (r = 0) to the contact line (r = 1) with 
So = 0.002, obtained by numerical solution of Eqs. (I42|l and 



height at r is /i(r, Z _1 (r)), and the volume concentration 
is exactly C max . Thus substituting the expression for the 
height profile (Eq. (pTTjO . we find 

S(r) = C max h(r, rHr)) = ^=(1 - r 2 )(l - ^(r)). 

(50) 

Alternatively, via the shock front function l(t) the density 
profile E can also be expressed as a function of time 

S(i) = ^p(l-/ 2 (t))(l-i), (51) 

representing the growth dynamics of the deposit. We 
shall use the same functional notation E for both the 
areal density profile E(r) and the growth dynamics E(t) 
interchangeably without confusion. 

A direct numerical solution for l(t) (Eq. (I42[l ) yields 
an areal density profile E (r) shown in Fig. 1121 and an 
areal density growth E (t) shown in Fig. [T3j The areal 
density profile shows that almost all the deposit occurs 
near the contact line with a pronounced peak followed by 
a steep decrease and then a long flat tail with negligible 
deposit. The areal density growth dynamics E (t), how- 
ever, shows that a significant amount of deposit occurs 
throughout the whole drying process with the maximal 
density much less pronounced in the time profile than 
the spatial profile. This is understandable from previous 
results (Fig. [TTj) as the shock front l(t) remains near the 
contact line for the most of the total drying time. 

We now study these profiles analytically via solving the 
solving for the shock front l(t). In the small So limit we 
seek the solution to Eq. (142|) via asymptotic expansion. 
Since < t < 1, we expand l(t) in terms of t and thus 
find 

l(t) = l-\Vht-^V5~ot 2 -^VS~ot 3 + 0(t i ). (52) 




0.2 0.4 0.6 0.8 1 
t 



FIG. 13: Areal density growth dynamics 2E(t)/(#iC ma x) with 
So = 0.002. The time domain is essentially from to 1 — So- 



We look for the position r m and the time t m — I (r m ) 
for the maximal areal density near the contact line. Nu- 
merical simulations show 

1 - r% = 0.31... V^, C = 0.67.... (53) 

Using first order approximation in t: l(t) « 1 — (l/4)yJot, 
we find 

1 - r W _ 0.125...V^o, t$ = 0.5. (54) 
Second and third order approximations in t yield: 

l-r£> =0.18...V^, t« = 0.55..., (55) 

1 - r (3) = 0.22...V^o, t$ = 0.59.... (56) 

It seems that different orders of approximation all show 
a contact line distance 1 — r m proportional to y/So but a 
maximal deposition time t m independent of <5o, and that 
with higher order of approximation the resulting 1 — r m 
and t m converge to the numerical results Eq. ([53")) . We 
thus conjecture that, in the small So limit, the distance 
between the maximal areal density and the contact line 
scales with y/fio, and the time at which the maximal areal 
density occurs is a constant portion of the total drying 
time. We shall confirm this conjecture later. 

We also note that since throughout the drying process 
the majority of the deposit occurs near the contact line, 
we can analyze the deposit spatial profile E (r) and the 
growth dynamics E (t) separately in different asymptotic 
regimes. 

B. Asymptotic Time-Independent Deposition 

We consider two cases where the time evolution of the 
deposition characteristics is negligible cither by focusing 
on the depositing process in the early time regime or by 
controlling the experimental conditions to assure time 
independence. 
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1. Deposition at Initial Drying Stage 

The initial drying stage corresponds to the asymptotic 
regime t/tf — » or equivalently tf — > oo. Similar to 
the previous study [2^ | , we consider a typical evaporative 
deposition scenario where there is a large pool of solution 
(tf — > oo), and we are only interested in the deposition 
process near the contact line (r — > 1) up to a time t/tf 
close to zero. 

In this approximation the time dependence of all the 
characteristic quantities is suppressed. Thus with t — > 
and r — ► 1 

Mr,*)«|(l-r- 2 )«^(l-r-), (57) 



J(r) 




t 



FIG. 14: Sketch of a possible realization of time-independent 
geometry. Solvent is injected into the drop at the center from 
below to balance the loss due to evaporation. 



v(r,t) u(l,0) = v , 



5*5 = 



Co 
c 

^ ' 1 1 1 n > 



and 



(58) 
(59) 

(60) 

(61) 
(62) 



The spatial profile (Eq. ([61])) is trivial as it is just the 
time-independent drop height profile (Eq. ([57]) ) multi- 
plied by the maximal concentration C max . The growth 
of deposit from the contact line is linear in time with the 
proportional coefficient determined by the initial condi- 
tions of the drying drop. 

Eq. (|6"Tj) also shows, at the initial drying stage, 



l(t) « 1 - y/S v Q t. 
The deposit profile satisfies: 

E(r) » 9iC maK (1 - r) , 



S(t) « ^C max (1 - Z(i)) « BiCn^y/SoVot. 



dE(r) 
dr 



(63) 



The slope of the deposit profile near the contact line only 
depends on the initial contact angle and the critical con- 
centration C max . We note that in this case the initial 
flow velocity Vq is finite at the contact line, and this is 
not a general property. We shall later discuss the case of 
the diffusion-controlled evaporation profile where vq di- 
verges, and we shall show the slope (Eq. ([63])) is actually 
independent of the evaporation profile. Intuitively, this 
property is not difficult to understand. Since the time- 
dependence is suppressed, the slope of the areal density 
at the contact line should be identical/proportional to 
the slope of the initial drop height profile when the dry- 
ing process begins, which is exactly the initial contact 
angle 0j. 



2. Deposition in Time- Independent Geometry 

An evaporating drop with time-independent geome- 
try can be achieved by replacing the evaporating fluid 
as illustrated in Fig. ([14"]) . Solvent is injected into the 
drop from below at the drop center with solvent mass 
J Q R 2i:rJ(r)dr per unit time so that the total amount of 
solvent in the drop is constant in time. The manipulation 
also allows minimal perturbations to the drop shape and 
the flow velocity field and thus achieves a drop height 
profile that is approximately time- independent. 

The hydrodynamics in this case is thus stationary, and 
with uniform evaporation profile 



h(r) = h(r,t = 0) = ^(l-r 2 ) 



v(r) = v(r, t = 0) = — r. 



(64) 



(65) 



Using the same method (Eqs. (03}, J34|), and ([38])) . we 
can show that the volume concentration profile in the 
transport region is: 

ft, ,n n exp (-t/2) - r 2 exp j-t) 

C(r,t) = C 5 , 66) 

1 — r z 

where Co is the initial uniform volume concentration. We 
show the volume concentration profiles C(r, t)/Co at four 
different times in Fig. [T5] As opposed to the time- 
dependent geometry, the concentration increases with 
time near the contact line but decreases close to the drop 
center, and there is a crossover. 

The spatial areal density profile is trivial: 



£(r) = C max / l (r) = 
and the deposit growth dynamics is 



E(t) 



' i '- max 



(i - i\t)) , 



(67) 



(68) 



where the evolution of the shock front l(t) is given by the 
pipe models Eq. ([5]). 
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FIG. 15: Volume concentration profile C(r,t)/Co (Eq. J66}) 
in the time-independent geometry at t = 0.2, 0.4, 0.6, and 
0.8, with 8 = 0.002. 



C. Density Profile in Spatial Asymptotic Regimes 

1. Near the Contact Line 

We now analyze the areal density profile in spatial 
asymptotic regimes. We first consider the limit l(t) — > 1. 
In addition we require, 



5 (i-*r i/2 «i, 



(69) 



that is, we consider in the time regime not too late in 
the whole drying process as we shall show later that the 
shock front reaches the drop center exactly at 1 — Sq. 
Evidently this is satisfied for any t < 1 for sufficiently 
small Sq. 

According to the previous results for the deposit pro- 
file E (r) and the growth process E (t) (Fig. [12] and Fig. 
[TB"]) . almost all the deposit occurs throughout almost all 
the drying process in this asymptotic regime. Thus this 
regime contains most characteristics of the deposition 
profile. In particular, we want to look at where (the 
contact line distance (1 — r m )) and when (time t m ) the 
maximal areal density occurs, and to show explicitly the 
scaling properties with respect to v^o- 

In this asymptotic regime Eq. ([42]) is reduced to 



dl 
dt 



(1-t) 



-1/2 



- 1 



4(1 -i) 



1-P 



(70) 



With the initial condition Eq. (143 p . l(t) can be solved 
for: 



(i-i(t)Y 



So 



2(l-t)~ 1/a + ln(l-t)-2l. (71) 



We substitute Eq. (p} into S(t) (Eq. ([51])) and look 
for the maximal value. We find the time t m satisfies the 
following equation 



21n(l-i ro ) + 3(l-i m ) 



-1/2 



Thus 



l(t, 



0.67, 



l-0.31..V*o 



and 



E (r TO ) = 0T1... 



(72) 



(73) 



(74) 



(75) 



The maximal areal density occurs roughly at 2/3 of 
the total drying time and is independent of the density 
ratio 5q when it is small. The distance from the contact 
line to the maximal areal density 1 — l m and the maximal 
density itself, on the other hand, scale with y/do. These 
results are consistent with our preliminary findings via 
asymptotic expansion. 

We can show the scaling properties more explicitly. 
Near the contact line r — ► 1 and 1 — r 2 = (1 + r)(l — r) w 
2(1 - r) (Eq. dH), and thus E(r) « ^C max (l - r)(l - 
Z _1 (r)). We introduce the scaled quantities 



E = 



and 



d=l-f 



1-r 



and Eq. (fBTj) shows 



1 - t 



E 

T 



(76) 



(77) 



(78) 



Substituting Eq. ([78]) into Eq. ([71]), we find 




(79) 



Eq. (|79p implicitly gives the scaled universal areal den- 
sity profile E as a function of the scaled contact line dis- 
tance d (Fig. \W\i . 

Eq. ([79"]) cannot be solved explicitly, and we instead 
look for the asymptotic power laws in different limits. 
We first consider the limit d — ► 0. On the right side of 
Eq. (|79p . the term containing logarithm dominates, and 
thus 



de 



4^ 



(80) 



The areal density increases with the contact line distance 
linearly. Eq. (|80p is consistent with our previous findings 
(Eq. dUD). 
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FIG. 16: Scaling property of the areal density profile for small 
So- The scaled areal density is defined as E = 2£/ ' (6iR\fSo); 
the scaled radial distance is defined as d = (1 — r) /y/So. The 
spectrum of curves going downward corresponding to respec- 
tively 5 = 0.1, 0.05, 0.01, 0.002, 0.001, and 0.0005. 




0.2 0.4 0.6 0.8 1 
t 



FIG. 17: Scaling property of the areal density growth dy- 
namics for small So- The scaled areal density is defined as 
£ = 2S/(6iRs/So). The spectrum of curves going down- 
ward corresponding to respectively So = 0.1, 0.05, 0.01, 0.002, 
0.001, and 0.0005. 



We next consider the tail of the profile in the limit 
d — > oo and E — > 0. In this limit (d/E) 1 / 2 is the dominant 
term, and thus 



E ~ -d~ 3 . 
4 



(81) 



In terms of the original varibles E and r, Eq. (|8ip can 
be rewritten as 



1 



•(l-r-r 



(82) 



The areal density profile is thus proportional to <5q along 
the tail as opposed to y^o toward the contact line. 

Similarly we can show the scaled areal density E as a 
function of time t explicitly: 



1 -t 



2(l-tr 1/2 + m(l-t) 



1/2 



(83) 



Eq. (fB"3"|) gives the universal dynamics of deposit growth 

(Kg. nzi. 



2. Toward the Center of the Drop 

Toward the center of the drop I — > 0, Eq. (|42[) takes 
the form 



dl 



I 6 (l- 1) 



-1/2 



dt 4(l-t)i_j (i_ t )-V2- 
We find in this case 



/ = A l-S (l-t) 



-1/21 1 / 2 



(84) 



(85) 



where A is a constant and can not be determined without 
proper boundary conditions. 

Thus toward the center of the drop l(t) oc l(t) (Eq. 
and as a result 



1(1 



5$)= I l-^)=0. 



(86) 



The shock front arrives at the center of the drop at time 
(1 — <5 2 )i/, and the whole horizontal transport process 
stops before the drop dries up. 
The areal density in this case is 



E(r) cx 6[ 



A 4 (l-r 2 ) 



(A 2 



(87) 



The density profile toward the drop center is proportional 
to <5q in contrast with the density near the contact line 
where it is proportional to v^o- This is consistent with 
our previous findings (Eq. I|82[)). For small So majority 
of the deposit thus occurs near the contact line. 

We also note that toward the center of the drop (I — > 
0) it is essentially in the very final stage of the whole 
drying process (t — * (1 — <5q)). Practically some model 
assumptions may not hold at this stage. For example, the 
drop height profile may no long decreases linearly in time 
However, the scaling property of E (r) with respect 
to 6% (Eq. I[87J)) may still be valid. 



V. DEPOSIT PROPERTIES WITH 
DIFFUSION-CONTROLLED EVAPORATION 

Evaporation profiles determine the hydrodynamics, 
such as the flow velocity field, and thus affect the deposit 
properties. Different evaporation profiles are achievable 
in experiments, and in some cases the resulting deposit 
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patterns change dramatically with evaporation condi- 
tions [31( . The uniform evaporation profile we have con- 
sidered is nonsingular while the diffusion-controlled evap- 
oration profile that has been generally assumed and stud- 
ied in the literature diverges at the contact line. We shall 
apply the pipe models to the deposition process with the 
diffusion-controlled evaporation in this Section to assess 
the effects of evaporation condition on general deposition 
properties. 



A. Hydrodynamics and Concentration Profile 

We derive the hydrodynamics and solute transport dy- 
namics of an evaporating thin drop with the same trunca- 
tion concentration C max but with the diffusion-controlled 
evaporation profile. 

Assuming vapor density above the liquid-vapor inter- 
face obeys the diffusion equation and taking into account 
the quasi-equilibrium shape of the thin drop, the evapo- 
ration flux defined as amount of mass per unit projected 
area per unit time is 



and 



J(r) oc (1 -r 2 )- 1 / 2 , 



(88) 



where we still use the dimensionless representation: r — > 
r/R and t -> t/t f with < r < 1 and < t < 1. 

The height profile takes the same form as with the 
uniform evaporation 



h(r,t) = ^(l-r 2 )(l-t), 



(89) 



though physically in the dimensional domain the im- 
plicit total drying time tf is defined with respect to the 
diffusion-controlled evaporation flux as 



p 27r/i(r, 0)rdr 
2nJ(r)rdr 



(90) 



Substituting Eqs. ([88]) and ([89]) into the condition of 
local solvent mass conservation (Eq. I]28p). we solve for 
the velocity field 



v(r,t) = -{l-t)- 1 - 
4 r 



(l-r 2 )- 1/2 -(l-r 2 ) 



(91) 



The height profile (Eq. 1(59])) and the velocity field (Eq. 
(|9"T]0 thus completely determine the hydrodynamics. 

To derive the solute concentration profile, we use the 
same method as with the uniform evaporation (Eqs. (|34[) 
and (|38p ). The previously defined initial position func- 
tion ro(r,t) (Eq. ([3"3"]) ) and the concentration C(r,t) 
in the transport regime have exact though complicated 
forms with the diffusion-controlled evaporation: 



r 2 ) 3/2 = l-(l-(l-r 2 ) 3/2 ) (l~i) 3/4 



(92) 



C(r,t) = C (l-t)- 1 / 4 (l-r 2 )- 1/2 



x 1- 



1- 1- r 



2\3/2' 



U-i) 3/4 } 



1/3 



(93) 



The pipe models can be applied to derive the equation 
of motion of the moving shock front 



di 



1-6' 



(94) 



where v(r,t) is given by Eq. ([9Tp . and S = C/C max with 
C(r,t) given by Eq. ||S1SJ|. _ _ 

At the contact line dl/dt — > — v(l, 0)v^o/(l — v^o)- It 
is finite for the uniform evaporation, but divergent in the 
diffusion-controlled case (Eq. ([91"]) ). We shall show later 
that this divergence does not affect the deposit properties 
near the contact line. 

As in the case of the uniform evaporation we can also 
introduce the virtual front l(t) by imposing the condition 
C{l(t), t) = C max , where C(r, t) is given by Eq. ([93]). The 
terminal time when the virtual front arrives at the center 
of the drop is given by the condition C(0, t) = C max , and 
we find 



t = 1 



o . 



(95) 



Thus the whole transport process stops at (1 — 5g) of 
the total drying time for the diffusion-controlled evapo- 
ration as opposed to (1 — Sq) for the uniform evaporation. 



B. Deposit Density Profile 

For the case of the diffusion-controlled evaporation, the 
areal density profile £(r) and the growth dynamics E(t) 
are completely solved by Eqs. ([50]), ([51]), and ([94]). We 
show the numerical results for the density profile and the 
growth dynamics in Fig. [T5] and Fig. [T5] The corre- 
sponding results in the case of the uniform evaporation 
are also shown for comparison. 

The diffusion-controlled evaporation thus yields more 
pronounced deposit density near the contact line as the 
singular flow field (Eq. ([9"T]) ) carries more solute particles 
toward the contact line at early time. For the areal den- 
sity profile, two asymptotic regimes can again be identi- 
fied as detailed later, with majority of the deposit occurs 
near the contact line. For the deposit growth dynamics, 
however, significant amount of deposit occurs throughout 
the whole drying process. Thus the two different evap- 
oration conditions lead to density profiles and growth 
dynamics with qualitatively similar properties. 

We now study the asymptotic regimes quantitatively. 
Toward the contact line, we consider the limit r — > 1 and 
6 — > So <C 1- The equation of motion for the shock front 
(Eq. ([M])) reduces to 
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and 
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FIG. 18: Deposition profiles with 5 = 0.002. The solid line 
corresponds to the case of the uniform evaporation. The dot- 
ted line corresponding to the case of the diffusion-controlled 
evaporation. 




FIG. 19: Deposition growth dynamics with So = 0.002. The 
solid line corresponding to the case of the uniform evapo- 
ration. The dotted line corresponding to the case of the 
diffusion-controlled evaporation. 



dt 



a- 



1/3 



-t)l- 



(96) 

With the condition 1(0) = 1, Eq. ([9"6"| admits a solution 
of the form: 

(1_;) 2 = ^ f 1 (l^xf^x^dx. (97) 

3 J(l_t)3/4 

Using the scaled density profile E (Eq. ([76")) ) and the 
scaled contact line distance d (Eq. (177)) '), we can rewrite 
Eq. ([97J as 



1 



(t/df 



(98) 



1 - t 



(l-t)3/ 4 



(1 



xf^x^dx. 



(99) 



Eqs. ([9"B")) and (JM)) thus give the universal density profile 
and growth dynamics in the limit <5o — > 0. 

Subjecting the growth dynamics Eq. ((99)) to the condi- 
tion of vanishing derivative with respect to time, we find 
that t m , when the maximal areal density occurs, satisfies 
the equation: 



2(1 -t m ) 
3 



(1 



x^x-^dx 



'(l-t m ) 3/4 

i(l-(l-< m ) 3/4 ) 1/3 (l-U 3/4 = 0. 



Eq. poop can be solved numerically, and thus 
t m = 0.50... 



(100) 



(101) 



It is worth noting that t m is not exactly 0.5, and Eq. 
(|100|) is not exactly symmetric about t m though the re- 
verse might seem true from Fig. 1191 

Substituting t m into the shock front solution Eq. ([57)) 
and the areal density profile Eq. ([50)) , we find the contact 
line distance for the maximal density 



1 — r m = 0.33... v<5( 
and the maximal areal density 

£ m = 0.17.. .V^C n 



(102) 



(103) 



We also look for the power law governing the tail de- 
cay of the scaled universal profile E. In the limit d — *■ oo, 
the right side of Eq. (|98[) must diverge and the integrand 
must be asymptotically dominated by x -4 / 3 in the lower 
limit (E/J) 3 / 4 -> 0. Thus Eq. ^ reduces to the follow- 
ing power law: 3 

£ ~ cT 7 , (104) 
and in terms of the original variables £ and r 

£ ~ 8iC„ 



°o 



(1-r)' 



(105) 



In comparison with the uniform evaporation profile, a 
more pronounced maximal areal density occurs earlier at 
about 1/2 of the total drying time as opposed to 2/3 of 
the total drying time, but at about the same distance 



3 Alternatively, we can derive the power law exponent ex- 
plicitly by using the formula f a (l - x) 1 /3 x -4/3 dx _ 

-3 (2(1 - p) 1 / 3 + P F?[2/3, 2/3, 5/3, p]) /(2p 1 /3). 
F^[2/3, 2/3, 5/3, p] is a hypergeometric function with 
[2/3, 2/3, 5/3, 0] = 1 and Ffp/3, 2/3, 5/3, 1] = T (1/3) T(3/5). 
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from the contact line. The areal density profile near the 
contact line still obeys the -^Jo-scaling law. Away from 
the contact line, however, the density profile decays much 
faster and is proportional to 8q as opposed to 6q for the 
uniform evaporation. 

We are for the first time able to derive the power-law 
decay (Eqs. l[52|) and (|105p ) of the deposition profile to- 
ward the center of the drop. Although we have based the 
derivation on specific deposition mechanism and simple 
conditions, these asymptotic results reflect the behavior 
of the shock front at a late time, and thus may reveal gen- 
eral properties of the evaporative deposition phenomena. 
We shall come back to this point later in Discussion. 



VI. COMPARISON WITH POPOV'S MODEL 

As mentioned in Introduction, Popov Q applied the 
same truncation criterion by Dupont to a similar problem 
with different conditions and assumptions. Unlike our 
analysis so far, Popov studied the solute transport pro- 
cess with diffusion-controlled evaporation and analyzed 
the final deposit pattern in terms of the thickness profile. 
These differences are not essential, however, as we shall 
discuss below. Most importantly Popov assumed that 
the deposition region interferes with the transport region 
and the hydrodynamics of the evaporating drop. In par- 
ticular, the formed deposit with finite thickness alters the 
boundary conditions under which the drop shape is de- 
termined. In Popov's model the drop shape is thus trun- 
cated by the deposit: the drop shape in the deposition 
region is identical to the thickness profile of the formed 
deposit while in the transport region it is described by a 
spherical cap elevated by the amount of deposit thickness 
at the boundary. As a result, the velocity field, deter- 
mined by conservation of the solvent mass, also depends 
on the deposit thickness, and the deposition process thus 
couples with the transport process. This assumption, 
though more realistic, makes the underlying mathemati- 
cal structure more complicated, and a complete solution 
is not available. Popov showed the asymptotic deposition 
properties near the contact line under these assumptions. 
We shall compare our findings with these results to shed 
some light on the model robustness. 



A. Thickness Profile 

We have described the final deposit pattern in terms of 
the areal density profile. To implement this description, 
we assume that in the deposition region, where horizontal 



4 It is consistent with the previous results for the terminal time 
(Eq. I I95H ) when the whole transport process stops. This can be 
shown by considering the concentration exactly at the center of 
the drop. 




FIG. 20: Sketch of the two descriptions. The shaded area 
is the deposition region. The dotted line outlines the height 
profile of the evaporating drop that decreases in time. The 
solid line shows the thickness profile of the formed deposit. 
A: description in terms of the areal density profile where so- 
lute particles move vertically in the deposition region and the 
height profile of the region follows that of the drop; B: de- 
scription in terms of the thickness profile where the thickness 
of the deposit is described by H(r). 



transport is stopped, solute particles can still move ver- 
tically so that the deposit thickness decreases along with 
the drop height profile and eventually reaches zero. At 
the same time, the volume concentration diverges and is 
replaced by the areal density profile. This zero-thickness 
assumption is not essential as long as the deposition re- 
gion does not interfere with the geometry and hence the 
hydrodynamics of the evaporating drop. Alternatively we 
can interpret the final deposit in terms of the thickness 
profile. In this description solute particles cannot move 
vertically either once the volume concentration reaches 
Cmax- Thus the thickness profile of the deposition region 
does not change with time once formed and the volume 
concentration is uniformly C max throughout the deposit. 
The evaporating drop, on the other hand, still maintains 
its spherical-cap equilibrium shape with height profile de- 
creasing in time, and the evaporation flux and the flow 
velocity field are not affected by the deposit. Essentially 
these two descriptions correspond to two different deposit 
dynamics in the deposit region during the drying pro- 
cess: in the areal density description, the deposit height 
profile is identical to the drop height profile h(r, t) and 
vanishes at the total drying time; in the thickness descrip- 
tion, the density profile is characterized by the function 
h(r, Z _1 (r)) and is independent of time. We sketch these 
two descriptions in Fig. [201 

The two mathematical formulations for the deposit 
profile are essentially the same as the solute mass across 
the thickness at each position in the deposition region 
is constant without horizontal transportation. The areal 
density profile is given by E (r) = C max h(r, l^ 1 (r)) (Eq. 
pi)) ), and thus 

ff(r) = £M =/l(r . r i (r)) (io6) 

(--'max 

is the corresponding thickness profile. Since H(r) is pro- 
portional to E(r) with the constant coefficient C max , our 
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previous findings for E (r) are readily applicable to H (r) 
with almost no modifications. 

The mathematical equivalence is based on the same un- 
derlying physical assumption that the deposition process 
does not interfere with the transport process and should 
be distinguished from Popov's model. This assumption 
may seem even farther from reality in the thickness de- 
scription as finite accumulation of the deposit almost 
certainly alters the velocity field and hence the trans- 
port process. However, as we shall show later deposition 
properties are actually independent of the specific phys- 
ical assumptions in some asymptotic regimes. 



B. Asymptotic Deposit Properties 

We first consider, in the limit t — ► 0, r — ► 1, the slope 
of the deposit thickness profile dH/dr. With £ given by 
Eq. O, 

^ = -r^c max (i - rV)) - ^ max (l ~ r 2 )^-^ 

= -fliCUc - ^iCUx(l - r 2 ) (j^j . (107) 

In the limit the second term on the right side of Eq. (|107|) 
is proportional to h(l,0)/v(l,0) at the contact line and 
vanishes. Thus near the contact line in the early drying 
stage dYi/dr — > — #iC max and thus dH/dr — > — 6i, which 
is independent of the evaporation profile (compared with 

Eq. HMD)- 

We next consider the growth dynamics of the deposit 
thickness H{t) = £ (t) /C max at the early drying stage. 
Since H (r) as shown above is proportional to (1 — r) in 
this limit, thus 

H{t)=H(r)^9i(l-l{t)). (108) 

For the uniform evaporation, it is given by Eq. (|62|) . For 
the diffusion-controlled evaporation we note that near the 
contact line dl/dt — * — \/Sov(l,t), and thus in this limit 

l-i(*)«V^(l-»-o(l,t)), (109) 

where ro (r,i) is the initial position function (Eq. (|92p ). 
With r = 1, t — * 0, and ro — * 1, the leading order term 
of H(r) (Eq. is 

if^VW^r, (HO) 
using (1 - rl) » 2 (1 - r ) and 1 - r » (1/2)(1 - (1 - 

4)3/4)2/3 „ ( 3t )2/3 /2 7/3_ 

The asymptotic form Eq. (|110p is exactly the same as 
the result derived in Popov's case Q- This sameness is 
not a surprise. In solving the system of coupled equa- 
tions, Popov argued that the solutions can be expanded 



in terms of a/^o- In the lowest order the equations de- 
couple, 5 and it physically corresponds to the condition 
that the deposition process and the hydrodynamics are 
separated. Thus our result Eq. (|110[) (as well as the 
assumption of decoupling) is the lowest order approxi- 
mation to Popov's findings near the contact line at the 
early drying stage. 

The later drying stage when t is large, however, was not 
well defined in Popov's results. Popov suggested that, in 
the limit 5q — > 0, the thickness profile H(r) monotoni- 
cally increases with t and eventually approaches a limit 
H m that is proportional to v^o toward the end of the dry- 
ing process (1 — t -C 1). The thickness profile thus ends 
abruptly with a "vertical wall" at its inner side. Although 
our finding of the maximal thickness H m = £ m /C max has 
the same v^-dependence (Eqs. (|75|) and (I103P ). the uni- 
versal properties of the profile show that H m is formed 
at about 2/3 (uniform evaporation) or 1/2 (diffusion- 
controlled evaporation) of the total drying time and the 
thickness decreases continuously afterward following a 
power-law decay (Eqs. ([82]) and (|105p ). Furthermore, 
our results clearly show two different asymptotic regimes 
and the thickness profile is instead proportional to 5q 
(uniform evaporation) or Sq (diffusion-controlled evapo- 
ration) in the late drying stage toward to center of the 
drop. Neither the v^-vmiverealhy of the thickness pro- 
file near the contact line nor the power-law decay toward 
the center of the drop was apparent in Popov's findings. 

The coupled system Popov studied is thus fully re- 
ducible to the simpler conditions we have assumed in the 
small 5q limit. This coupling between the deposition and 
the hydrodynamics is further resolved near the contact 
line in the limit r — > 1 (this is consistent with the v^o - 
scaling property of (1 — r)). It is thus fairly reasonable to 
argue that the evaporative deposition phenomenon is ro- 
bust against specific assumptions and conditions in this 
asymptotic regime. Additionally as we have shown that 
the temporal properties of the process are actually Jo- 
independent, by assuming the decoupling we are able to 
specify a universal time scale when the maximal deposit 
thickness forms and to extend the truncated deposition 
model to a separate new asymptotic regime toward the 
center of the drop. 

Furthermore, it is worth pointing out that even 
Popov's treatment in [7J] only partly addresses the cou- 
pling between the formed deposit and hydrodynamics, 
and a few important dynamical aspects are still miss- 
ing. First, the deposit-dependent drop shape alters the 
boundary conditions under which the evaporation flux is 
determined and hence further changes the velocity field. 
Second, the velocity field in the deposition region should 



5 The solution Eq. JTlOj corresponds to the zeroth-order in \/3o 
in the asymptotic expansion Q. The explicit dependence on 
^/So in Eq. IllOt . however, appears via a different mathematical 
route. In Dupont model, it is due to the boundary condition Eq. 
J 16H at the contact line. 
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be certainly affected by the deposit. If the deposit takes 
up a nonzero fraction of the volume, a given evapora- 
tion flux will take the same amount of available solvent 
per unit time but with a larger effective volume in the 
deposition region, and a larger velocity is required to re- 
plenish this solvent. However, our findings suggest that 
by controlling the initial concentration ration Sq all these 
couplings might be effectively suppressed near the con- 
tact line at the early drying stage and may not alter 
the deposition properties. Also, further considerations of 
these couplings are important for experimentally testing 
the model. We shall discuss this issue later. 



VII. COMMENTS ON DEPINNING AND 
MULTIPLE-RING PATTERNS 

In the evaporative deposition model, we have identi- 
fied one key deposition mechanism: horizontal transport 
of solute particles truncated by maximal allowed volume 
concentration. This truncated dynamics is robust in the 
sense that it admits a solution with a moving shock front 
that has a definite functional form independent of the dy- 
namical characteristics of the problem. The shock front 
velocity only depends on local values of the flow veloc- 
ity, the pipe geometry, and the volume concentration in 
the transport region. Furthermore, the truncated trans- 
port mechanism leads to a characteristic deposit profile 
near the contact line that is independent of evaporation 
conditions and dynamical couplings in the early drying 
stage. 

Although practically evaporative deposition phenom- 
ena are rich and there must be many complex underlying 
mechanisms, we now discuss how the truncated transport 
dynamics, as a single deposition mechanism, might yield 
general deposit profiles under some typical and control- 
lable conditions. 



A. Depinning Configuration 

For an actual evaporative deposition process Q , in the 
late drying stage the contact line of a evaporating drop 
can depin and retreat toward the drop center before it 
is pinned again, and the deposition process is restored 
at the newly-formed contact line. The slip-stick process 
may occur many times before the drop tries up and thus 
yield a multiple-ring pattern [25|, [H, [H, [HJ ■ 

There is still no consensus on a general depinning 
mechanism, which may not exist after all. A conven- 
tional explanation 0, @, l26l ca lls for the Young-Dupre 
equation at the contact line j35|: 7;„ cos# + 7 s ; — j va = 0, 
where 9 is the contact angle, and ji v , 7 S /, and j va are the 
surface tensions at liquid-vapor, solid-liquid, and vapor- 
solid interfaces. The horizontal component of the surface 
tension along the liquid-vapor interface at the contact 
line 7;„ cos 9 acts as the depinning force that increases 
as the drop progressively flattens toward the late drying 



stage. j s i also contributes to the depinning force, but 
it primarily depends on the local deposition and is less 
affected by the geometry. 

Admittedly this explanation can be problematic in 
that it predicts an extremely narrow range of deposit 
thickness and drop configurations when the maximal hor- 
izontal depinning force is achieved and the contact line 
depins. In reality, however, the depinning occurs over a 
wide range of configurations. To account for this smooth 
change, the conventional explanation thus must be com- 
plemented by other mechanisms that are still lacking. 
We shall in the following consider within the framework 
of the conventional theory and cautiously keep its limi- 
tations in mind. 

The depinning phenomenon thus clearly involves the 
coupling between the deposition process and the hydro- 
dynamics of the evaporating drop: the drop geometry 
and therefore the contact angle 9 must be determined 
from the boundary conditions defined by the formed de- 
posit pattern to account for the change of the depinning 
force (Fig. 1.2A). To continue our discussion of depin- 
ning within the model established in the previous Sec- 
tions, we must adopt the view that our results based on 
the assumption of decoupling can be applied to the cou- 
pling system in the lowest order of So, and we also must 
interpret the deposit pattern in terms of the thickness 
profile H. 

Deegan showed experimentally [f| that the depinning 
time td is about 0.4 to 0.8 (in the unit that the total dry- 
ing time tt = 1) and it depends on the initial volume con- 
centration <5o (from 0.001 to 0.01). We have found that in 
the small Sq limit, the maximal areal density and hence 
the maximal deposit thickness occurs at time t m ss 0.67 
(Eq. ([75]) ) for the uniform evaporation and 0.50 (Eq. 
(jlOip ) for the diffusion-controlled evaporation. Our t m 
thus falls in the range of the depinning time found in ex- 
periments. Geometrically the maximal deposit thickness 
occurs when the contact angle 9 at the rim of the deposit 
becomes zero and the depinning force 7;„ cos 9 achieves 
its maximal value. 

To derive t m formally, we note that H(t) = 
£(<)/C max = (9 l /2)h(l(t),t) (Eq. (EU)), and thus the 
condition dH/dt = demands 



dh dh dl 
dt dr dt 



(111) 



Eq. (jllip does not depend on model assumptions as 
h(r,t) can be the drop height profile under any condi- 
tions. In particular, h(r, t) can be the solution to the drop 
geometry coupled with the deposition process, dl/dt, on 
the other hand, must be given by the pipe model re- 
sults Eq. ([5]). Eq. (|111[) gives the general condition and 
scheme to solve for t m when the maximal thickness occur 
within the truncated transport model. 

However, t m is only an approximation for the depin- 
ning time td- The depinning may occur once the contact 
angle is small enough but not necessarily zero, and as a 
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result td < t m . 6 More importantly, our finding of t m , 
as opposed to Deegan's results, does not depend on the 
initial concentration (Fig. [17]) when So is small. This de- 
pendence may be included in the higher order solution in 
So when the coupling between deposition and geometry is 
considered. The dependence can also reflect the depen- 
dence of j s i on So- Nevertheless our finding of t m within 
the truncated transport model is consistent with empiri- 
cal evidence and may be instructive for further studies. 

Also, for the complicated depinning process the above 
explanation based on Deegan's experimental results is 
not universal. Other mechanisms, such as the tension 
within the formed deposit [36[ , the kinetic processes con- 
trolled by the viscous stress in the drop, or the diffusion 
relaxation in air, may affect the depinning under certain 
conditions and thus lead to different interpretation of the 
depinning time. 



B. Formation of Multiple- Ring Patterns 

The depinning process can lead to multiple-ring de- 
posit patterns: the retreating contact line is rested at a 
nearer position toward the drop center and the deposi- 
tion process continues at this newly formed contact line. 
In Refs. [1H and (26[, a slip-stick theory was established, 
where the Young-Dupre condition was incorporated into 
the equation of motion for the contact line. As a result, 
the contact line was found to follow a sinusoidal motion, 
and this led to a very dense multiple-ring pattern. The 
authors in Refs. [25[ and [33j considered controllable 
multiple-ring formation in nonconventional geometries. 
Some of these empirical evidence may need explanations 
that are different from Deegan's framework [a]. 

Following our discussion on the depinning process, we 
consider an ideal model for the multiple-ring formation. 
For the diffusion-controlled evaporation, we identify the 
depinning time td = t m «l/2 (2/3 for the uniform evap- 
oration), i.e., the contact line is released when the con- 
tact angle becomes zero. The width of the deposit ring 
is identified as (1 — l m ) (Eq. (74]) ). The thickness is rep- 
resented by the maximal thickness H m = S m /C max (Eq. 

The newly pinned position R' can be determined by 
considering the change of the solvent volume. Initially 
for small contact angle di, the total volume of the drop 



V 



-R 9; 



(112) 



We assume the contact line moves fast enough to the new 
position so that the time it takes can be ignored compare 



As noted by Popov , the contact angle may become negative 
when the drop is actually overhanging against the deposit rim. 
In this case > tf, and the depinning process must be due to 
different mechanisms. 



to the time when it is pinned. Since the evaporation flux 
is stationary, the volume of the drop when the contact 
line is pinned again is thus V(l — td) ~ V/2 (V/3 for the 
uniform evaporation). For slow evaporation, the newly 
pinned drop must still maintain a spherical cap geometry. 
The contact angle at the new contact line must still be 
approximately 6i since it is determined by the Young- 
Dupre condition. Thus we find 



R' 



1/3 



R^0.8R, 
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and R' w 0.7R for the uniform evaporation. 

We next consider the width and the thickness of the 
new deposit ring. The deposition process via truncated 
transport is restored at the newly formed contact line, 
and the deposit growth there must follow what we have 
found at the initial contact line. Since R(l — r m ) cx Ry/So 
(Eq. (TD) and H m cx Ry/% (Eq. (75])) (we now have 
added R to the dimensionless expressions), for the new 
ring we replace R with R' and replace So with concentra- 
tion ratio at the new contact line: S' w C(R', td)/C max w 
2.2<5o in the case of the diffusion-controlled evaporation 
profile (Eq. (38])) (2.4o"o for the uniform evaporation). 
Thus the width R'(l — r' m ) and the thickness H' m of the 
new ring satisfies: 



R'(l - r' m ) « l.LR(l - r m ), 



H' 



1.1H„ 



(114) 
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For uniform evaporation, they are 1.2iZ(l — r m ) and 
1.2H m respectively. The thickness and height of the sec- 
ondary ring are thus comparable to those of the initial 
ring. 

There can be further depinning events in this ideal 
model. Following the above reasoning, for the n-th order 
ring (the initial ring is the zeroth order with the radius 
R) toward the drop center, the time it starts to form 
(the depinning time of the (n — l)-th ring), the radius, 
and the concentration ratio are 



.(n-l) 



1 71— X 



k=0 



1 
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and 
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C o2n/3 1 

^max * ' — 1 

Thus the width and the thickness of the n-th ring are 
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FIG. 21: Photographs of the deposit patterns in Deegan's 
experiment JjJ. The left column shows the entire deposit for 
initial volumn concentration (a) 0.01, (b) 0.0025, (c) 0.0013, 
(d) 0.00063. The right column shows the close-ups. (Courtesy 
Robert Deegan) 



n/2 



/27n/6 
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The results for the uniform evaporation can be derived 
similarly. 

Empirically deposit patterns with major multiple rings 
were observed by Deegan Q. In Fig. for two initial 
concentrations (a) 0.01 and (d) 0.00063 there are major 
higher-order rings up to the third order with width and 
thickness comparable to the zeroth order ring at the ini- 
tial contact line. Our above analysis in the highly ideal 
case is at least consistent with the observed pictures. For 
the other two initial concentrations (cases (b) and (c)), 
however, major higher-order rings are not observed. 

Between the major rings Deegan's results clearly show 



some fine complex patterns. In our ideal case, these fine 
patterns are formed by the solute particles in the region 
where the depinned contact line sweeps through between 
two pinning events. Thus they must be formed within 
a much smaller time scale compared to the major ring 
formation according to our assumptions. 7 Formation of 
those fine structures, as a result, may not be explained 
within the framework of the truncated transport dynam- 
ics considered in this work. Instead, other mechanisms, 
such as the dynamics and instabilities of the dewetting of 
thin films [38, 39, 40] where capillary forces are important 
and long-range interactions between solute particles are 
relevant [U |42j , are needed to account for the observed 
patterns in those regimes. 



VIII. DISCUSSION 

A. Experimental Aspects of the Model 

Experimental tests involve two aspects of the model. 
First, Dupont model is based on some assumptions and 
conditions that require specific experimental design and 
realization. Second, testability depends on the major 
predictions made by the model. We shall discuss these 
in turn. 

The decoupling between the deposition process and the 
transport process is the major assumption. This ideal 
case can be approximated under certain experimental 
conditions. 

One of the key parameters in Dupont model is the al- 
lowed maximal concentration C max . C max should depend 
on the size of the solute particles, interactions between 
the solute, the solvent, and the substrate. C max could be 
also affected by other dynamic factors, such as the flow 
velocity. Empirical data for C max is limited, but naively 



can be deduced from the initial volume concentra- 



tion by measuring the volume of the major deposit rings 
and comparing it with the initial volume of the entire 
drop. 

The assumption of decoupling requires C max be small. 
Any solid particles must become immobilized if their con- 
centration is high enough. For simply shaped particles, 
the immobilization occurs at volume fractions of order 
unity. However, the immobilizing volume fraction can be 
much smaller than unity if the solute consists of tenuous 
objects like polymers, especially those that form a phys- 
ical gel (3f| . Many organic solutes such as food may well 
be of this sort, and seem most appropriate for the model. 

Another aspect of the coupling is that the formed de- 
posit may alter the drop geometry which in turn changes 
the evaporation flux and its generated flow field. One way 



7 We estimated in [37| that the typical velocity for a free standing 
dewetting film is about the order of 1 m/s. 
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to bypass this obstacle is to realize the uniform evapora- 
tion profile as sketched in [2^ | . The uniform evaporation, 
unlike the diffusion-controlled case, is maintained by the 
peripheral environment but not the drop itself. The cou- 
pling due to change of the evaporation condition is thus 
effectively avoided. 

Furthermore, some kinetic effects such as Marangoni 
flows and viscous stress may deviate from the Dupont 
model assumptions and the fundamental lubrication ap- 
proximation. Convective Marangoni flows can be avoided 
by maintaining uniform experimental conditions such as 
the temperature and by reducing the thermal effect due 
to evaporation. Viscous stress effect can be also reduced 
by slowing the evaporation. Thus by making the envi- 
ronment nearly saturated with the solvent vapor, these 
kinetic effects can be made arbitrarily small. 

Our work allows one to compare the predictions of the 
Dupont model with experiment quantitatively in several 
new ways. Our main prediction is about the normal den- 
sity profile. There are three aspects that might be tested: 
the leading edge, the time and position of the maximal 
deposition, and the trailing edge. 

The leading edge prediction states that the areal den- 
sity should increase linearly with distance from the con- 
tact line, and that the slope is proportional to the square 
root of the initial concentration ratio. This prediction is 
identical to Popov's results and is therefore not a par- 
ticular test of Dupont model. The maximal deposition 
time, on the other hand, is a strong test and comes out 
very clearly from our findings. 

The power-law decay of the areal density at the trailing 
edge is also a strong test. On the qualitative level, this 
decay appears consistent with some experiments but not 
others. Often with colloidal solutes, such as in some of 
Deegan's experiments, one sees a sharp inner edge with 
no apparent decay at all. This is not a surprise. As 
noted above, Dupont model is not completely applica- 
ble with compact solid solutes. Moreover, a continuum 
model such as Dupont 's should not work completely when 
the solute particles have a size comparable to the drop 
thickness h.Foi non-colloidal solutes, there should be a 
decay that is consistent with our predictions for C ma x 
small enough. 

For quantitative prediction, one should measure the 
density or thickness of the deposit as a function of dis- 
tance as well as drying time. The function of distance 
can be straightforwardly measured for the dried deposit. 
The function of drying time can be measured by monitor- 
ing and recording the whole drying process. The shock 
front velocity might be also measured this way. 



B. Generalizations of the Model 

We have not considered any kinetic effects such as vi- 
sous stress and diffusion in the work. These effects will 
modify the physical assumptions and mathematical for- 
mulations. For example, we implement a step function 



to model the deposition condition with respect to the 
threshold concentration C max . Horizontal diffusion and 
forces between various media may smooth the transition 
to immobility. Mathematically this effect can be mod- 
elled by replacing the step function with a weakly con- 
vergent series, and the truncated transport equation thus 
becomes nonsingular. This is an interesting mathemat- 
ical problem. Physically we want to point out that the 
smooth effect is not as necessary as it might seem. When 
Sq is small, the concentration in the transport region rises 
so fast toward the contact line that the transitional region 
with partial mobility would be very small. 

We already suggested that the power-law decay (Eqs. 
([52"| and (|f 05[) ) of the trailing edge could be independent 
of specific assumptions and derivation, and may actually 
reveal general properties of the evaporative deposition 
without depinning. To be more specific, as already partly 
discussed in the previous work [29( asymptotic properties, 
such as the power-law decay, might be robust against a 
broad range of specific dynamic processes, but are in- 
stead dependent on the singularities that characterize the 
dynamic conditions. In the case of the evaporative depo- 
sition, there are two sources of singularities encountered: 
the vanishing of the drop height and the divergence of 
the diffusion-controlled evaporation profile. Both occur 
at the contact line. The exponent of -7 in the power- 
law decay of the density profile might be uniquely de- 
termined by these two singularities. This thought may 
prompt further studies of the robustness of the power-law 
against the drop geometry and other singular conditions. 
For instance, an evaporating drop pinned in an angular 
region studied in [29( has a singular geometry profile at 
the contact line dependent on the contact line distance 
from the vertex of the angle as well as the opening angle 
of the region. An exact relation between the decay ex- 
ponent and the characteristic singularities would give a 
decaying profile that continuously varies along the con- 
tact line and is controllable by the opening angle. More 
detailed asymptotic analysis is needed toward this goal. 

Although in this work the simple geometry (spheri- 
cal cap) and the regular evaporation profiles (uniform or 
diffusion-controlled) yield consistent and rich deposition 
phenomena, controllable nonconventional geometry and 
evaporation (and so is the flow velocity field) are impor- 
tant for creating complex deposit patterns [13L l33l . |34| . 
Rigorous theoretical studies on these patterns are chal- 
lenging. For example as pointed out by Popov Q, de- 
termining the evaporation profile with nonconventional 
boundary geometry is one of the major obstacles for rig- 
orously solving the deposition problem. Furthermore, 
some drop configurations and velocity field, such as the 
negative contact angle, might lead to the formation of 
higher-order deposit rings even before depinning hap- 
pens. 

The pipe models we have established in this work cap- 
ture the basic characteristics of the truncated transport 
dynamics. It is robust in the sense that it admits a mov- 
ing shock front, of which the equation motion only de- 
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pends on local characteristic quantities in a definite func- 
tional form. Besides the evaporative deposition prob- 
lem, pipe models may be relevant to other systems with 
similar properties, such as the jamming in granular flow 
[43l l44j . Further studies are needed to adapt the pipe 
models to such systems. 

In addition, some aspects of the pipe models may in- 
vite further attention. For example, to understand the 
behaviors of the shock front at singularities of a pipe and 
to characterize and utilize such behaviors require our fur- 
ther efforts. We have only considered in this work the sin- 
gularity associated with the vanishing height of the pipe 
(Eq. (|16p ). and in the application to the evaporative de- 
position problem this singularity corresponds to the con- 
tact line and amounts to the condition of self-consistency 
(see Appendix for details). However, the formulation Eq. 
(fT6f is far more general than the application. Unlike the 
case of evaporative thin drops where the height profile 
and the flow velocity field are coupled due to the mass 
conservation of the solvent, In our derivation of the pipe 
models, the height profile and the flow velocity field are 
allowed to have independent dynamics. Thus the singu- 
larities can rise either dynamically as those characteristic 
quantities evolve or statically in a controllable way. How 
to describe the shock front with dynamical singularities 
is an interesting but challenging problem. 



IX. CONCLUSION 

We have studied in this work a highly idealized evap- 
orative deposition model in which different dynamics is 
fully separatable. In addition, the model is scalable in 
the sense that the length scale, defined in terms of the 
initial drop radius, and the time scale, defined in terms 
of the total drying time, are irrelevant in model charac- 
teristics. Yet this model has yielded rich and solvable 
deposition properties that are consistent with empirical 
evidence. 

In retrospection, we have so far fairly successfully ad- 
dressed the two goals suggested in the Introduction Sec- 
tion. 

The first goal is to break the problem down and to 
prioritize different components. We divide the contact 
line deposition problem into two parts: phenomenology 
and mechanism. The former includes various conditions 
such as evaporation profiles, capillary-dominant equilib- 
rium drop shape, and flow velocity field. For the lat- 
ter part, among all the complexities we have identified 
and focused on the truncated transport dynamics as the 
unique deposition mechanism. Our findings suggest that 
this unique mechanism is effective and robust for a fairly 
wide range of phenomenologies to produce consistent de- 
posit patterns. 

Our analysis has revealed a few asymptotic regimes 
with different characteristic scales in space and time. 
These regimes are largely determined by the phe- 
nomenologies but not the mechanism. The truncated 



transport dynamics is applicable to some of the regimes, 
and asymptotic deposit properties and scaling laws are 
identified. But some other important regimes may only 
be addressable by other mechanisms that go beyond the 
scope of this work. 

The second goal is to find the appropriate mathemati- 
cal structure, which is largely determined by the mecha- 
nism. We have found the "pipe models" and the charac- 
teristic moving shock front solution useful structures to 
describe the truncated transport dynamics and the de- 
position problem. It has a universal functional form and 
depends only on local conditions specified by the phe- 
nomenologies. It is adaptable and generalizable. 

In summary we thus have found: transport dynam- 
ics truncated by the allowed maximal volume density, 
combined with the hydrodynamics under lubrication ap- 
proximation of a thin drop evaporating regularly with 
uniform or diffusion-controlled evaporation flux, gives a 
simple, tractable, and yet rich enough picture for the gen- 
eral contact line deposition phenomena. A few asymp- 
totic regimes both in drying time and radial distance are 
identified and analyzed. We are able to predict in cer- 
tain cases deposit profile and growth properties that are 
consistent with empirical evidence and previous findings. 
These properties are largely robust against physical con- 
ditions and theoretical assumptions, and are instructive 
to further applications. 
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APPENDIX A: DIRECT DERIVATION OF THE 
SHOCK FRONT VELOCITY 

The shock front l(t) moves from the contact line (r — 
1) toward the center of the drop (r = 0). Behind the 
boundary 1 < r < l(t) is the deposition region; in front 
the boundary, which is in the transport region, solute 
particles are carried by the flow toward the boundary 

(Fig.ua). 

We consider the change in total solute mass in the 
region l(t) < r < 1. Initially the amount of solute mass 
is 

M = [ 2nC h{r, 0)rdr. (Al) 
Ji{t) 
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C{r,t) 



h(r,0) 
.h(r,t) 



r = \ 



l(t) 



r=0 



and solve for dl/dt: 



dl 



dm/dt 



dt 2nl [C h{r, 0) - C max h(r, t)] - dm/dl 
I 



4(1-*) 



C m „(l-I 2 ) 

^[(l-t)- 1 / 2 -; 2 ] 



(A5) 



FIG. 22: Sketch of the definition of Mo and Aft. 



Eq. 



51) can also be written in the form: 



For a position r with l(t) < r < 1 when the moving 
front reaches r at time Z -1 (r), the volume concentration 
at r becomes C max and the height at r is h(r, Z — 1 (r)). 
The amount of solute mass near r at time Z _1 (r) is 
thus 2irrC max h(r,l~ 1 (r))dr. Since there is no horizon- 
tal movement of solute particles in the deposition region, 
this amount of mass does not change in time. Therefore 
at time t the total amount of solute mass in the region 
l(t) < r < 1 is 



2TrC mliK h(r,r 1 (r))rdr. 



(A2) 



Hi.) 



Conservation of solute mass demands that the change 
in total mass M t — Mq must be equal to the amount mass 
transported through l(t) from time zero to time t, i.e., 

-l ,i 

27rC max /i(r, r^r^rdr- / 27rC ft(r, 0)rdr = m(l{t). 
H(t) Ji(t) 

(A3) 

m(l(t), t) is defined by Eq. (|3"4"|) . We differentiate against 
t on both sides of Eq. fAlij) . 

- 2,lC max h(r,t)- + 2«lC h(r,0)- = ^ + ^ 

(A4) 



dl 
dt 



-v(l,t) 



1-5' 



and we thus have re-covered Eq. (f42|) . 



(A6) 



The initial condition Eq. (|4"3"1) can be proved via 
asymptotic expansion. Near the contact line t is small, 
and we can write l{t) = 1 — at + 0(t 2 ). Substituting this 
expression for l(t) on both sides of Eq. |42|) and retaining 
only the lowest order of t in the limit t — > 0, we find 



We solve for a: 



<5 (l+4a) 



a 



*), 



(A7) 



(A8) 



Eqs. (|A5[) and (|A8|) completely define the equation of 
the motion for the shock front. 



t=o 
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The initial condition Eq. (I43|) is thus proved 
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